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Abstract 

This work provides a complete description of entanglement properties between 
electrons inside coupled quantum molecules, nanoestructures which consist of 
two quantum dots. Each electron can tunnel between the two quantum dots 
inside the molecule, being also coupled by Coulomb interaction. First, it is 
shown that Bell states act as a natural basis for the description of this physical 
system, defining the characteristics of the energy spectrum and the eigenstates. 
Then, the entanglement properties of the eigenstates are discussed, shedding 
light on the roles of each physical parameters on experimental setup. Finally, a 
detailed analysis of the dynamics shows the path to generate states with a high 
degree of entanglement, as well as physical conditions associated with coherent 
oscillations between separable and Bell states. 

Keywords: Entanglement, quantum molecules, quantum information. 

PACS: 68.65.Hb, 03.67.-a, 03.67.Bg 


1. Introduction 

In the last two decades, the interest on the physical behavior of semicon¬ 
ductor materials has increased due to their potential applications in quantum 
computing BE]- In particular, semiconductor quantum dots have been proved 
to be an ideal candidate for the codification of quantum information. In these 
systems, a qubit can be defined using the charge |3] and the spin of the particle 
confined mini, as well as excitonic states Bi- An interesting system is 
a quantum molecule (QM) which is a nanoestructure consisting on two quan¬ 
tum dots separated by a barrier that can be tunneled by charged particles, as 
electrons mmmm 

In 2003, Hayashi et al. demonstrated the coherent manipulation of a charge 
qubit in a QM m- The quantum nanoestruture is built by using metal gates 
for confining charges on a bidimensional electron gas in GaAs/AlGaAs. The 
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same group implemented a quantum device with two QMs coupled electrostat¬ 
ically but isolated by conduction m- In this system, the resonant tunneling 
current through each molecule is influenced by the charge state of the sec¬ 
ond molecule m- More recently, researchers on the Quantum Device Lab 
(ETH Zurich) and Universite of Sherbrooke (Quebec, Canada) demonstrated 
successfully the coupling between the same type of device with a microwave 
resonator m- The experimental data (together with a theoretical treatment) 
shows that the coupling between a quantum molecule and the resonator follows 
a Jaynes-Cummings type interaction. In another interesting work, quantum os¬ 
cillations between three quantum states in a Si/SiGe double quantum dots are 
controlled by high-speed voltage pulses DU- Measurements of the transconduc¬ 
tance permitted the observation of this oscillations between different coupled 
quantum states induced by two different pulse profiles. Both works show the 
potential applications of this kind of heterostructures on quantum information 
processing, opening the possibility of the definition and control of quantum bits. 

From the theoretical point of view, Fujisawa et al. [18] demonstrated that, 
in the experimental setup mentioned above, it is possible to implement some of 
the two-qubit gates, including the Bell gate at a very specific choice of physical 
parameters. Other theoretical works have explored aspects related to quantum 
correlations and decoherence [19, [20:, also considering restricted physical condi¬ 
tions. Due the importance of the role of entanglement in quantum computation 
protocols, a complete analysis of the entanglement properties of this particular 
system, with potential application of quantum information processing, is crucial. 

Here, it is developed a careful exploration of the entanglement properties 
of electrons in coupled QMs, considering the realistic experimental setup of 
Ref. [14], treated as a closed system. The results show that Bell states are the 
key behind the physical behavior of the electrons inside the coupled QMs. Also, 
it is demonstrated how the careful control of physical parameters as tunneling, 
electronic energy offsets (detunings), and Coulomb coupling, can be used to cre¬ 
ate highly entangled states. First, using analytical calculations together with 
numerical simulations, the characteristics of the energy spectrum and the eigen¬ 
states of the Hamiltonian are explored. The results show that the eigenstates 
correspond to Bell states under specific physical conditions. Second, the effect of 
physical parameters is mapped, focusing on the generation of highly entangled 
states. The “beats” on the dynamics of entanglement are verified and explained 
as a competition between two different frequencies associated with the physical 
couplings. Physical requirements to obtain coherent oscillations between sep¬ 
arable states and a specific Bell state are discussed. All results are obtained 
considering realistic values of physical parameters. 

This paper is structured as follows. Section [2] contains the description of 
the theoretical model together with the definition of the measurement of en¬ 
tanglement, concurrence, used in this work. Section [3] is devoted to explore the 
characteristics of the energy spectrum and the eigenstates from the point of view 
of entanglement properties. Section [4] is reserved to the discussion of dynamics, 
focussing in the obtention of Bell states using temporal evolution, considering 
the effects of tunneling and detuning between electronic states inside each QM. 
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The work is summarized in Section [5] 


2. Model 


The analysis performed here is based on an actual experimental setup with 
two coupled QMs mm, each with two quantum dots separated by a potential 
barrier which can be tunneled by electrons mmmm- The quantum dots 
for each QM are arranged horizontally, being coupled to a source and a drain 
of electrons, used to charge and discharge the QM. Coulomb blockade and the 
physical design of the nanoestrutures guarantee that there is only one extra 
electron and, effectively, one electronic level for each quantum dot. Each QM 
can be treated theoretically as a two-level system: the electron can occupy 
the left (right) dot | L) (|i?)). Even if the two QMs are designed in order to 
inhibit intra-molecules tunneling, both subsystems are still coupled because of 
the Coulomb interaction between the electrons on each molecule. From the 
point of view of quantum information, this correspond to the implementation of 
two qubits, one on each QM, which are coupled in such a way that the physical 
system is bipartite. 

This system is well modeled by the Hamiltonian m written as 

H = \ E ( e ‘d° + Ai<r«) + J -o W (1) 


where Si parameters describe the energy offsets of the electronic levels, also 
known as detunings, and A i are tunneling rates inside each QM. Last term is the 
Coulomb interaction between electrons from different QM, which is quantified 
by J. Considering a\ = \L) u (L\ - \R) U (i?| and ct* = | L) u (i?| + \R) U (L\, the 
matrix form of Hamiltonian 0, written in the positional basis of the coupled 
molecules (|MiM 2 ) = | LL) , | LR ), | RL ), \RR}}, is given by: 


H = 






( 2 ) 


where e s = s i + £ 2 and Sd = e\ — e 2 . As measurement of entanglement, we use 
the concurrence, an efficient measurement for bipartite 2 0 2 systems HQ [22], 
which is defined as 


C(p) = max {0, Ai - A 2 - A 3 - A 4 } . (3) 

Here A^ are the square roots of the eigenvalues, in decreasing order, of the 
matrix R = pp. The p operator is obtained through the spin-flip operation 
of the complex conjugate of the density matrix of the 2 0 2 system so p = 
(cfy(S)cry)p^(o-y(S)0'y). If C(p) = 0, the state of the bipartite system is separable. 
On the other hand, when C(p) = 1, the state is a maximally entangled state. In 
the following sections, the theoretical model and the entanglement measurement 
are used to explore the effects of physical parameters on both, the entanglement 
properties of eigenstates, and the entanglement dynamics. 
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3. Entanglement properties of eigenvectors 

Consider the Bell states given by [23] 

l*±> = ±={\RL)±\LR)), 

l*±> = ±={\RR)±\LL)). (4) 

These are eigenvectors of the Hamiltonian © when Ai = A 2 = 0 and 61 = 62 = 
0 and also form a complete basis for 4 O 4 Hilbert space. Using the ordering 
Ob = (|T_) 5 13>-), |4/ + ) , the Hamiltonian © becomes: 
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where A± = (Ai db A 2 )/2. Written in the Bell basis, it becomes evident that 
tunneling and Coulomb couplings split the Hilbert space in two 2 x 2 sub¬ 
spaces, described by matrices M±, with the matrix D coupling both of them. 
Rewriting Hamiltonian ([l]) using the Bell basis is advantageous to explore the 
rich entanglement behavior of the system: it separates the effect of tunnel¬ 
ing A_(A+), which is responsible for the oscillations within 2 -fold subspaces 
|T_) (|T + )) and |4>_) (|4> + )), from the effect of detunings and responsible 
for inter-subspaces mixing. This new approach is ideal to study the entangle¬ 
ment properties of eigenstates using both, analytical and numerical calculations. 


3.1. Analytical solution at e 1 = e 2 = 0 

An analytical solution is obtained by considering the resonant electronic 
levels so ei = e 2 =0. Notice that, for this condition, the 2 x 2 matrix D in 
Eq.(j5| is null so the four Hamiltonian eigenvalues are given by 


E^=T 


J 2 + 16A| 


with the corresponding eigenstates written as 


ip**} = ey (At) +ri- T i$=f>) 
^ (|$t> + I*=f>) • 


Here, the new quantity <4 T T is defined as 


at =t 


j ~ 




16A| 


4A 3 


(6) 


( 7 ) 


(8) 
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while T = l/yjl + U 2 is the normalization term. The eigenstates are now a 
superposition of two elements on Bell basis which depends on the tunneling 
coupling. It is straightforward to check the validity of this solution by con¬ 
sidering the limiting case of tunneling suppression (A x = A 2 = 0). For this 
condition, the function (j = =p) goes asymptotically to zero when A^ 0 
(i = =f) and, consequently, T —>► 1 so the set given by {) , 1^+) , \^P-) 5 1^+)} 
coincides with {|T_), |<f>_), |T + ) , |<F + )}. 

The analytical solutions, Eqs.( |6|7| ), are useful to explore the physical conse¬ 
quences of fixing equal tunneling rates (A 2 = A x = A). Under this condition, 
the matrix M_ in Hamiltonian © becomes diagonal, while M + remains non¬ 
diagonal with A + = A. Depending on the relations between parameters e\ and 
£ 2 , the system falls into one of the following behaviors: 

• At full resonance, when £ 2 = £1 = 0, the eigenstates correspond to Eq.([7|), 
with two of them corresponding to the Bell states |T_) and |<h_) with 
energies — J/4 and J/4 respectively. The others are given by |V^) with 
energies 

• Considering £2 = £ 1 , which means £d = 0, the eigenstate 1^-) remains as 
|T_), a maximum entangled state, while couples with subspace “+” 
by the action of detuning e 3 . 

• If £2 = — £ 1 , then = 0, the eigenstate becomes the Bell state |4>_) 
while the other three eigenstates are coupled. 

3.2. Numerical results 

In this subsection, the entanglement properties of the coupled QMs are ex¬ 
plored in a wider scenario: first, the action of detunings is studied when the 
system has equal tunneling rates, A_ = 0. Then, the effects of both tunneling 
terms, A + and A_, are quantified. 

3.2.1. Equal tunneling rates 

Fig. [T] shows the numerical results for concurrence of the eigenstates of 
Hamiltonian © as a function of detunings £$, considering equal tunneling rates 
(A_ = 0 and A + = A). In this calculations, two different choices of tunneling 
parameters are used: A = J/16 (weak tunneling), Figs. [lja-d), and A = J/4 
(strong tunneling), Figs, [lje-h). The Coulomb coupling is fixed at J = 25 /ieV, 
as reported in experiments PUTS]. For all cases, the eigenstates are shown in 
ascending order of energy so {10) , 11) , 12) , |3)} are ground state and the first, 
second and third excited states respectively. 

First, the analytical results for the three types of conditions discussed above 
are confirmed with the numerical calculations: at full resonance, eigenstate 
|0) and 13) correspond to 1^-) an d 1^+) respectively while |1) = |T_) and 
12) = |4>_) are maximally entangled states, fact confirmed by checking the 
coefficients of the eigenstates calculated by the numerical simulation and shown 
by the bright yellow line over the condition £2 = £1 in Figs. 0b) and (f). For 
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the condition Sd = 0 , | 1 ) remains in |T_), while at £2 = —£ 1 , the second excited 
state 1 2 ) = |<F_), conditions associated with yellow lines seen in Figs. [ljc) and 
( g )- 

Out of the three resonant cases, the value of concurrence has an abrupt 
decrease due to the enhancement of coupling between subspaces associated with 
M±. The eigenstates |0) and |3) show a lower entanglement degree, if compared 
with | 1 ) and 1 2 ), with some differences that can be seen in Fig.[lja) and Fig.[lje), 
for state |0), and Fig. Bd) and Fig. 0 h), for state |3): at the full resonant 
condition and weak tunneling rate, |0) remains as \f>-) while |3) = both 

being superpositions with high concurrence value given by C(p) — 0.9. Then, 
following the condition £2 = £1 in Fig. [lja) and £2 = —S\ in Fig.[ljd), the value 
of concurrence decreases being C(p) ~ 0.5 at the two points corresponding to 
\£i\ = J/ 2. The decreasing value of concurrence is the signature of coupling 
between subspace M + with one of the Bell states: for |0), the coupled states are 
I'tp-) and |4>_) while in state |3) the coupling is between |?/+) and |T_). It is 
expected that the three choices for £i which favored highly entangled states as 
eigenstates of the system affect the behavior of the entanglement dynamics of 
this system. 

For values of £i out of resonance conditions, the four Bell states are cou¬ 
pled by both, the tunneling and detunings, and the mixing of basis elements 
explains the lower value of concurrence. Nevertheless, by comparing the two 
cases of tunneling, the values of concurrence increase when A increases, fact ex¬ 
plained because tunneling favored the coupling between internal states on each 
subspaces, which forms a superposition less sensitive to the action of detuning. 
Still, the numerical calculations show that the high degree of entanglement of 
superpositions with C ~ 0,97 obtained for A = J/16 at full resonance 

condition, Figs. [lja,d), decreases as tunneling rate increases, being C ^0,70 
for the same resonant condition in Figs. 0 e,h). 

It is interesting to note that, concerning detunings, there are two different set 
of parameters: the first set is defined by condition £i E [—J/2, J/2] where sub¬ 
spaces M + and M_ are coupled with each other. For the second set of values, 
it is verify that the detunings have the effect of exchanging the value of the 
energies associated with | 1 ) and | 2 ) for even and odd resonances, still coupling 
the four Bell states out of any resonant conditions. For higher values of A, 
Figs, [lje-h), the increase of value of tunneling rate has the effect of enhancing 
the coupling inside the subspace associated with M+, which causes the increase 
of the “area” with high values of concurrence, 0.45 < C < 1 as can be easily 
check by comparing (a) and (e) in the figure. Thus, the set of conditions for 
entangled states for A = J/16 expands covering a wide region, although the 
expansion happens around the fixed critical points £i = =t J/ 2 . 


3.2.2. Different tunneling rates 

Fig .[ 2 ] shows the behavior of concurrence for Ai 7 ^ A 2 , considering Ai = J/4 
and A 2 = J/ 8 . In this case, the eigenstates of the system are not Bell states 
at any condition of resonance. Inside the region with £i E [—J/2, J/2], it is 
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e 2 /J e 2 /J 



Figure 1: Concurrence of the eigenstates of Hamiltonian (TT1) as a function of £\ and £2 with 
J = 25 jueV. For Ai = A 2 = A = J/16: (a) 10), (b)|l), (c)|2) and (d) |3). For Ai = A 2 = 
A = J/4: (e) |0>, (f)|l>, (g) |2> and (h) |3). 
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Figure 2: Concurrence of the eigenstates of Hamiltonian (TTl) 
sidering Ai = J /4 and A 2 = J /8 for J = 25 fieV: (a) |0), (b) 


as a function of £\ and £2 con- 
|1), (c) |2) and (d) |3>. 


observed the effect of the coupling between subspaces M + and M_: around 
the conditions 82 = £1 and 82 = —£ 1 , there is a region with higher entangle¬ 
ment degree for all eigenstates, although the concurrence is higher for the states 
associated with M_, |1) and |2). This is explained because the effective tunnel¬ 
ing A_ is weaker than A+. Out of these limits, there is a region with a high 
degree of entanglement for eigenstates |1) and |2), a new feature if compared 
with Fig. [l] That means different tunneling rates expand the range of physical 
conditions for obtaining entangled states with a high degree of entanglement 
(C > 0.7), although Bell states are not longer attained. 

4. Dynamical generation of entangled states 

A path for creation of highly entangled states is the use of dynamical evo¬ 
lution. First requirement is to prepare the system in an initial state. At ex¬ 
periments, it is possible to manipulate the energy offsets to set a specific con¬ 
figuration on couple QMs. The process is highly controlled m and |T(0)) can 
be fixed as one of the states {| LL ), | LR ), | RL ), | RR)} of the positional basis. 
Tunneling parameters are defined by construction of the nanoestructure, also by 
additional gates that inhibit or enhance coherent tunneling. The setup is open 
to manipulation of the detuning parameters £*. The readout of the state after 
the evolution is performed manipulating the bias on QMs and measuring the 
contribution of one electron to the current M- Inside this context, the effects 
of the physical couplings on quantum dynamics can be studied using analytical 
and numerical simulations. 
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Figure 3: Concurrence dynamics considering l'F(O)) a | RL) with J = 25 ^ueV. (a) Concurrence 
as a function of time and Ai / J considering Ai = A 2 for 0 < t < 3 ns; (b-e) Concurrence 
evolution for 0 < t < 1 ns considering specific choices for A/J ratio. 


4-1. Tunneling effect on entanglement dynamics 

The first goal is to analyze the effect of tunneling rates, once the results 
obtained in Section [3] show the relevance of those coupling in the characteristics 
of the energy spectrum. A numerical approach provides insights about the 
dynamical preparation of highly entangled states. In Fig. |3ja), the obtained 
numerical results of concurrence are shown as a function of time and Ai/J 
ratio considering |T(0)) = | RL), J = 25 /ieV, Ai = A 2 and resonance condition 
given by 62 = e\ = 0. It is verified that the tunneling coupling has a strong 
effect on entanglement dynamics: the evolved state goes from separable (dark) 
to a highly entangled state (bright) after the evolution time t e which depends 
on Ai/J. This can be verified by comparing the four profiles in Figs. Ib-e), 
noticing that the time when the concurrence has a local maximum at times 
which depend on the value of Ai/J, also shown in Figs. [3jb-e). It is interesting 
to observe that the oscillations are not fully periodical. This kind of behavior 
is a signature of the dynamical competition between different couplings, each 
with a characteristic frequency. This is better seen in Figs. IF -e), where the 
patterns resemble the “beats” obtained with the superposition of two harmonic 
oscillators. 

To shed more light on the physical origins behind this concurrence beats, 
the time-dependent Schrodinger equation is solved considering the full resonant 
condition, 51 = £2 = 0. This can be done using the previous results for eigen- 
ergies, Eq.([6|, and eigenstates, Eq.([7|). If |T(0)) = \RL), it is straightforward 
to obtain the evolved state |T(£)) and then the populations for the states of 
computational basis, Pij = | (ij|T(t)) | 2 . The last are written as: 


Prl = 


1 

4 



sin sin Q_t 


2 
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( 9 ) 

( 10 ) 


with the new physical parameters Q± defined as 

p± \[j ~ 2 + 16A ± 
0± = T = -4- 


( 11 ) 


These two frequencies are behind the concurrence beats in Fig. [4] Because 
Q± are functions of the tunneling parameters A±, the mechanism behind the 
rich entanglement dynamics shown in Fig. [3] is the coupling between Bell states 
inside the subspaces M±. That means the action of the tunneling can be used 
to manipulate the obtention of highly entangled states. 

There is an exact condition for coherent periodic oscillations between sepa¬ 
rable and Bell states: notice that states | LL) and \RR) are not populated for 
evolution times corresponding to 


nir 



( 12 ) 


with n being any integer number. At these specific times, the expressions Plr 
and Prl on Eq.(10) reduce to 

Prl = « jsin 2 fl_t e + [(-1)" + cosfl_£ e ] 2 j , 

Prl = | jsin 2 f \-t e + [(-1) 


) n -cos Cl-t, 


] 2 }' 


(13) 


In this way, the behavior of the function cos(^ + t e ) defines the condition for 
generation of Bell states |4/±) in Eq.Q. If its value is zero, the population of 
both states goes to 1/2, indicating that the system goes to a Bell state. The 
evolution time is also given by 


t e 


rrm 

2fil’ 


m odd. 


(14) 
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Figure 4: Populations and Concurrence dynamics for initial state given by l'P(O)) = \RL). 
Here J = 25 yeV, £2 = s\ and Ai = A2 being Ai/J = \/3/4 which fulfills condition given 
by Eq.( |15| > considering n = m = 1. (a) Populations of the four positional states \M 1 M 2 ) as a 
function of time; (b) Periodical oscillation of Concurrence. 


Comparing Eq.(12) and Eq.(14), both evolution times are equal if condition 
m = * s satisfied- In term of coupling parameters, the condition is equivalent 
to 

T <15) 

which has a real solution for m < 2 n. The populations and concurrence dy¬ 
namics are shown in Fig. [4j for physical parameters which follow the conditions 
given by Eq.(15). Notice that, at , the initial state | RL) exchanges population 
with \LR), to return again after two periods of oscillation. At times defined 
by Eq.([l4]), at times when the evolved state is a separable one, Bell states 
with concurrence value C = 1 are obtained as a superposition of the | RL) and 
| LR) states. Numerical and analytical results obtained by considering the other 
three states of positional basis \M 1 M 2 }, not shown here, establish that the ini¬ 
tial states given by | LR) or | RL) evolve to Bell states |T±), while Bell states 
|3>±) are dynamically created by setting | LL) and \RR) as initial states. 


4.2. Effects of detuning on concurrence dynamics 

The analysis of the concurrence behavior of eigenstates is a guide to check 
the effects of detunings given that, additionally to the full resonant condition, 
highly entangled states becomes eigenstates of the system at conditions given 
by £2 = £1 and £2 = —£i- When thinking about dynamics, the presence of such 
eigenstates at these specific conditions could favor the dynamical generation of 
Bell states. 

To obtain a full insight of the effects of detuning, numerical simulations of 
twelve different situations, considering as initial states the four elements \M\M 2 ) 
of positional basis and the conditions £2 = £1 and £2 = — £\ for £\ E [—J, J], 
chosen based on the analysis of spectrum as discussed in Sec. [3] The results are 
shown in Fig. [5] Each of the twelve physical situations fall in one of the three 
different patterns as follow: 
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• Fig. Ufa) : 

1. |T(0)) = | RL) and | LR) with £2 = £ 1 ; 

2. |T(0)) = | LL) and \RR) with £2 = —e 1 

• Fig.gb): 

1. |T(0)) = |i?L) and 52 = £ 1 ; 

2 . |T( 0 )) = |LL) and 62 = —£ 1 . 

• Fig.^c): 

1. |T(0)) = | LR) with £2 = — £ 1 ; 

2. |T(0)) = | RR) with £2 = £ 1 . 

Based on these results, we conclude that choices of detunings out of the full 
resonant condition do not break the link between the sets {| RL) , | LR}} and 
{| LL ), each with a specific detuning condition associated with oscilla¬ 

tions between a product of separable states and the corresponding Bell states. 
This is closely related to the characteristics of the eigenstates discussed in Sec. [3} 
the formation of Bell state |T_) (|<h_)), a superposition of | RL) and | LR) (| RR) 
and | LL)), is favored by the condition £2 = £1 (£2 = — £ 1 ). 

Nevertheless, the full resonant condition optimizes the coherent oscillations 
for all cases, as shown by the yellow regions between — J/4 < £\ < J/4 in 
Fig. [5|a). Out of this choice of parameters, the system falls into one of the 
asymmetrical patterns shown in Fig. §b) and Fig. [5jc) . Depending on the 
initial state, the creation of high entangled states is favored by negative or 
positive values of One feature common in the three patterns is the change 
of the behavior of the oscillations at values of \ei\ > J/2. This is also explained 
by the characteristic of eigenstates, which shows a different behavior for values 
out of £i G [—J/2, J/2], as show in Fig. [l] and Fig. [ 2 ] 

1 

0.8 
0.6 
0.4 
0.2 
0 

0 1 2 30 1 2 30 1 2 3 

t (ns) 



Figure 5: Three patterns for concurrence dependence of detunings Si and time, for tunneling 
couplings (Ai — A 2 ) and Coulomb coupling J so Ai/J = a/ 3/4, (a) |4^(0)) = (| RL) , | LR)} 
for £2 = £1 and |^(0)) = {| LL) , | RR)} for £2 = —£i\ (b) 1^(0)) = | RL) considering £2 = £1 
and |T(0)) = |LL) for condition £2 = —£ 1 ', (c) 1^(0)) = | LR) considering £2 = £1 and 
|T(0)> = \RR) for £2 = —£i- 
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5. Conclusions and perspectives 

Here, it is presented a careful analysis of entanglement properties of two 
electrons inside coupled QMs, focusing on the characteristics of the eigenstates 
and the quantum dynamics. The concurrence is used as entanglement mea¬ 
surement, once there are two coupled qubits defined in this particular system. 
First, it is demonstrated that the physical system is intrinsically linked to Bell 
(maximally entangled) states. Although the primary coupling connecting the 
QMs is the Coulomb coupling ( J) between electrons, the entanglement proper¬ 
ties of the eigenstates depends strongly of the tunneling coupling (AQ and the 
energy offsets (e*), also called as detunings . Specifically, the tunneling splits 
the 4x4 Hilbert space in two 2x2 subspaces of Bell states, while the role of 
detunings is to mix the two subspaces. The physical conditions for obtaining 
highly entangled eigenstates are fully described by analytical and numerical cal¬ 
culations, which includes specific parameter choices for obtaining Bell states as 
eigenstates. 

Second, it is explored the dynamics of the system in order to create highly 
entangled states. The results provide a full characterization of the role of tun¬ 
neling and detuning. When ressonant electronic levels (si = 0) are considered, 
it is shown that the tunneling coupling is the key behind the controlled gen¬ 
eration of Bell states. The apparition of “beats” in the dynamics is caused 
by competition between two different frequencies, which are functions of both, 
Coulomb and tunneling couplings. This beats can become coherent oscillations 
between positional states and a Bell state for a specific choice of ratio between 
tunneling and Coulomb couplings. Out of the resonant condition, this coherent 
oscillations are very sensitive to changes on the values of detuning. 

As future works, we want to search for a proposal of experimental measure¬ 
ment of the degree of entanglement in this system. In a theoretical proposal 
Emary [20 demonstrated that the process of decay of the qubits on coupled 
quantum molecules, in the long-time limit, induces rotations on each qubit 
as well as the direct decay. Based on this fact, the author propose a corre¬ 
lator which can be used as measurement of entanglement through a CHSH 
inequality [24]. The correlator is a function of concurrence as defined in this 
work. One of our perspectives for future work is to explore the definition of a 
different entanglement witness, looking for signatures of entangled states in the 
electronic current, generated by the system when pumped and drained. Another 
future research is the understanding the effects of dephasing and decoherence 
mechanisms over general dynamic, as well as the search for robust states. 
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